Superclimb of Dislocations and the Anomalous Isochoric Compressibility of Solid "^He 
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In the experiment on superfluid transport in solid '*He [PRL 100, 235301 (2008)], Ray and Hallock 
observed an anomalously large isochoric compressibility: the supersolid samples demonstrated a 
significant and apparently spatially uniform response of density and pressure to chemical potential, 
apphed locally through Vycor "electrodes". We propose that the effect is due to superclimb : edge 
dislocations can climb because of mass transport along superfluid cores. We corroborate the scenario 
by ab initio simulations of an edge dislocation in solid *He at T = 0.5K. We argue that at low 
temperature the effect must be suppressed due to a crossover to the smooth dislocation. 
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At present, the experimental search for supersolidity 
(proposed theoretically in Refs. [l|) in ^He focuses mostly 
on torsional oscillator experiments 0, and on attempts 
to detect pressure driven non-plastic flow So far, di- 
rect superflow through solid ^He has been observed only 
in the experiment by Ray and Hallock 0, Q . Although 
it was argued that the presence of a few liquid channels 
is compatible with the observations 0, 01, the absence 
of flow at temperatures above T 0.6K is a strong ar- 
gument against this scenario and in favor of superfluid- 
ity along dislocation cores or grain boundaries. Theo- 
retically, a superfluid dislocation network can manifest 
itself as a genuine superfluid or be in the Shevchenko 
state , characterized by anomalously low viscosity due 
to phase slips. In practice, the Shevchenko state might 
mimic superfluidity even at relatively high temperatures, 
T > O.IK, well above the actual transition determined by 
the dislocation density. 

The "UMass sandwich" setup of Refs. [1, Q is dif- 
ferent from the pressure driven cells 3. Superfluid 
''He is fed into the crystal through Vycor "electrodes", 
meaning that the chemical potential ^ is the physical 
quantity relevant to the external perturbation applied 
to the crystal. An insulating (i.e., non-supersolid) crys- 
talline groundstate has to be isochorically incompress- 
ible: X = {dn/dfi)Y = 0; that is. the density n of the 
crystal should demonstrate no response to infinitesimal, 
quasi-static changes of /i [§]. Indeed, as long as the cre- 
ation of single vacancies and interstitials is forbidden by 
finite energy gaps, the only way the density of a crystal 
can react dynamically to a small change in the chemi- 
cal potential, 6fi, is by creating/removing crystalline lay- 
ers. This requires nucleation times exponentially large 
in Thus, at temperatures much smaller than the 

vacancy /interstitial gaps, the isochoric compressibility x 
associated with thermally excited vacancies and intersti- 
tials is exponentially small. Consistent with these ar- 
guments, all non-supersolid samples of Refs. [1, Q have 
X — 0- two pressure gauges monitoring the solid showed 
no response to a change in /i by the Vycor electrodes. 



Supersolids have no vacancy (interstitial) gap [T^l and 
are thus genuinely isochorically compressible: X 7^ 0- 
One might argue that x should scale linearly with the 
superfiuid fraction ps since both are due to zero-point 
vacancies (interstitials). Given the extremely low value 
of Ps < 10~^ following from estimates based on the ob- 
served supercritical fiux (see Refs. for more details), 
one does not expect a noticeable x- However, a den- 
sity/pressure response to p. by several orders of magni- 
tude larger than expected is precisely what was observed, 
to which we refer as the effect of anomalous isochoric 
compressibility. Remarkably, the response was appar- 
ently spatially homogeneous, since two pressure gauges 
attached to two ends of the solid typically showed equal 
variations (but different absolute values; most samples 
were characterized by a static pressure gradient) [1, 

In this Letter, we argue that the microscopic phe- 
nomenon behind the effect of anomalous isochoric com- 
pressibility in the experiments by Ray and Hallock is the 
superclimb of superfluid edge dislocations, that is, climb 
controlled by superfluid flow along the core. Our idea is 
that significant and spatially uniform mass accumulation 
in the bulk of supersolid "^He is due to the synergy be- 
tween: (i) the presence of a superfluid network capable 
of delivering ^He atoms from Vycor electrodes to distant 
bulk regions and (ii) the presence of edge dislocations, 
whose superclimb is responsible for the density/pressure 
change. 

We corroborate our scenario by ab initio simulations 
which show that edge dislocation with Burgers vector 
along the hep C-axis has superfiuid core (we previously 
reported the superfluidity in the core of a screw dislo- 
cation [HI), and that it can climb in response to varia- 
tions of p. We argue that at low temperature the climb 
must be suppressed due to a crossover from a rough to a 
smooth dislocation [l^ . This prediction is a manifesta- 
tion of the structural evolution of dislocations with tem- 
perature, and is important for experimental validation of 
the scenario. While superflow is a necessary condition for 
superclimb, the dislocation must also have a finite den- 
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sity of jogs to allow for thrcshold-less climb. Otherwise, 
a finite gap A for creating dislocation jogs will protect 
the dislocation from shifting significantly in response to 
small variations in fi. 

The effect of anomalous isochoric compressibility is one 
of the novel properties emerging in the "quantum met- 
allurgy" [l^l context. These properties have long been 
discussed in the past; for example, it was speculated that 
quantum dislocations should be characterized by "thick" 
(roughened) cores due to zero-point motion [l^ . An im- 
portant role of quantum roughening of dislocations in the 
torsional oscillator response has also been proposed in 
Refs. [isl, U^. SupercHmb is a quantum analog of classi- 
cal high-T climb due to thermally activated flux of vacan- 
cies toward, away or along the cores (pipe diffusion) [l7| 
which adds (removes) atoms to (from) the extra plane 
forming the edge dislocation, so that the dislocation core 
shifts along the extra-plane direction. Obviously, at low 
T, the activated mass flow is exponentially suppressed 
and quickly becomes negligible. 

Apart from climb, dislocations can also glide. In 
Ref. it was shown that gliding dislocations (gliding 
does not require mass influx) are smooth at T = be- 
cause Coulomb-type interactions between shape fluctua- 
tions [l^, [l^ induce an energy gap AgUdo with respect 
to creating a pair of kinks in Peierls potential. Hence, 
threshold-less glide of a dislocation can effectively occur 
only at T comparable with AgUdo. This gap is also re- 
lated to shear modulus stiffening at low T [20[ . Similarly, 
dislocations have a gap A for creating a pair of jogs at 
T = 0, which leads to a suppression of climb (and x) at 
low T. The values of A can be quite different from Agiidc 
because the jog-anti-jog dcconfincmcnt couples to fluctu- 
ations of the superfluid density leading to an additional 
mechanism for the gap formation. 

Model of climbing superfluid dislocation. We introduce 
a coarse-grained description of an edge dislocation with 
superfluid core oriented along the X-axis in terms of the 
core displacement ?/(x,t) along the Y-axis (in the climb- 
ing direction) , perpendicular to the Burgers vector which 
is along the Z-axis. We proceed under the assumptions 
of small gradients and large displacements compared to 
the lattice spacing. Then, a coarse-grained density vari- 
ation 6n{x,t) translates directly into a coarse-grained 
variations Sy{x,t) oc Sn{x,t). The proportionality co- 
efficient is purely geometrical: adding one atom to the 
edge results in its displacement by a lattice period in the 
climb direction Sy{x, t) — a! and also in a density change 
Sn{x,t) = 1/a, where a is the length of the unit cell along 
the core. Thus, 6n{x,t) = £_Sy{x,t) with = 1/aa'. This 
relation implies that for a superfluid dislocation the core 
displacement Sy is the conjugate variable to the super- 
fluid phase if. The combined coarse-grained, low-energy 
effective action for superfluid and displacement degrees of 
freedom in the imaginary time description reads {h = 1) 



S = 



where the purely dislocation part of the action, Sd, is 
taken in the form of the Granato-Liicke string subject to 
Peierls potential [ll,[2l|: 
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with ni being the linear mass density of the core, Vd 
standing for speed of sound along the string determined 
by shear modulus G: « G/ni, and u denoting the 
strength of Peierls potential. In Eq.([21), the kinetic en- 
ergy oc y'^ is neglected in the low energy limit under the 
consideration. Full quantum mechanical description of 
the system based on calculating the partition function 
J DyDipeyii>{—S) will be presented elsewhere. 

Apart from the Peierls term oc u (not to be confused 
with the sine-Gordon term where the argument would be 
oc y{x')dx'), the quantized action ([I])-© is a standard 
harmonic (1 -I- l)-dimensional action. A renormalization- 
group analysis, similar to the one given in Ref. [l^ . shows 
that, at T = 0, the Peierls term has scaling dimension 
dim[u] = 2 regardless of the parameters of the system, 
even if the long-range deformation potential forces are 
ignored. This means that the Peierls barrier is relevant 
at T = and always leads to a finite gap A for the 
climb motion, i.e. the dislocation in its groundstatc is 
smooth. In such a state, the cosine term can be ex- 
panded in powers of y around some equilibrium position 
2/m = ma', m = 0,±1,±2,.... Accordingly, in the low- 
energy limit — when the gradient in the action ^ can 
be ignored — the action ^ reduces to the standard ID 
superfluid action (2^ 
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with g = u(27r)^/a'^. This action describes superfluidity 
with speed of sound vi = y/pTg/S, oc y/pTu and also a 
finite climb in response to Sfi: Sy ~ ^5fi/g. 

With increasing T, thermally excited jogs and kinks 
render Peierls potential less and less relevant, so that 
eventually it can be ignored. In this limit, the disloca- 
tion becomes rough, that is, similar to a free string [2l| . 
and the spatial gradient in Eq. ([2]) should be taken into 
account. The effective action ([T]) then becomes 



So = 



dr / dx 



dr / dx [-t^yip + {ps/2){d,^f - fi^y] + Sd,{l) 



Eq. ^ predicts an extremely strong quasi-static climb re- 
sponse: d^dy oc — (S/i determined by the length of a free 
dislocation segment L (the cross-linking distance in the 
network), so that a typical displacement Sy oc L^S^. This 
implies that the resulting specific compressibility is inde- 
pendent of the dislocation density « ^/L'^, provided the 
network is uniform over the whole sample. Indeed, the 
added amount of atoms per each "elementary" cube of 
the side L is aLSy oc L^Sfi. Thus, the added fraction 
of atoms per unit volume is independent of L. 



(4) 
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The superfluid component also demonstrates an un- 
usual behavior. The equation of motion for small oscil- 
lations reads 



(5) 



meaning that the spectrum of superfluid excitations is 
not sound-like anymore. It is described by a quadratic 
dispersion uj = ^Jr\q^ ^ where q is the momentum along 
the dislocation line. Full quantum mechanical descrip- 
tion of the crossover from the regime ([3]) to ([1]), (O in 
line with the approach of Ref. [l^] will be presented else- 
where. Here we point out two qualitative predictions 
of the model IT])-©: (i) suppression of the climb at 
T < A, and (ii) dramatic softening of superfluid phonons 
at T > A. 

Numerical results. Our ab initio Monte Carlo (MC) 
simulations were based on the worm algorithm (23j . The 
most important numerical finding of the present study is 
that edge dislocations with Burgers along the hep axis 
have superfluid cores in solid ''He. Our example is based 
on the dislocation with the core along the X-axis (and 
Burgers vector along the Z-axis) . Since the hep structure 
has two atoms in the unit cell, two extra half-planes are 
involved. Figures [T][H show snapshots of atomic positions 
in a typical MC configuration, along C-axis and along 
the core. Particles outside the circle, Fig. [21 were pinned 
to their classical lattice positions and provided boundary 
conditions for the simulation cell. The studied disloca- 
tion splits into two partials with the fee fault forming in 
between 
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The splitting is so large that it does not 
fit the simulation cell since one of the partials has moved 
all the way to the cell boundary. A direct simulation of 
the fee fault yielded an unmeasurably small (within our 
accuracy) fault energy < O.lK/atom, meaning that the 
splitting (proportional to the inverse of the fault energy) 
is indeed expected to be as large as > 150 — 300A . Cor- 
respondingly, physical properties of the both partials are 
essentially independent from each other. 

Under these circumstances we performed extensive 
simulations of a single partial attached to the fault. The 
rectangular simulation cell contained from 600 to 3400 
particles with periodic boundary conditions along the 
core. In perpendicular directions a boundary of pinned 
^He atoms surrounding a cylinder of radius R provided 
the necessary boundary conditions for the simulated sam- 
ple of solid ''He containing the partial dislocation at the 
center and the fault extending in the positive Y-direction. 
Depending on R, the number of actually simulated parti- 
cles varied from 270 to 1700. Superfluid properties were 
detected by observing winding exchange cycles along the 
cylinder axis (X-axis). The core response to changes in 
fi has been studied by tracing the position of the maxi- 
mum Y(fj,) of the columnar superfluid density map in the 
(Y,Z) plane [H. 

The core position exhibited strong continuous response 
to variations of fj,. The slope dY/dfi was larger in bigger 
cells indicating that at the simulated temperature T ~ 
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FIG. 1: (Color online) Columnar view of a typical MC con- 
figuration along the C-axis: filled red dots show atomic posi- 
tions; open blue circles indicate an ideal lattice; vertical solid 
green lines mark positions of the partial cores. The fee fault 
is between these two lines. The superfluidity occurs along the 
green lines. 
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FIG. 2: (Color online) Columnar view of the same MC con- 
figuration (as in Fig. [TJ along the cylinder axis. The circle 
marks the simulation cell where particle positions have been 
updated. The superclimb occurs in the "horizontal" plane. 
Green stars mark positions of the two partials. 



0.5 K it is controlled by the image forces provided by 
the boundary conditions. At fixed jj,, the configuration- 
to-configuration fluctuations of the core position were as 
large as several unit cells. Remarkably, the exchange- 
cycle map does not show any visible modulation with the 
lattice period in the Y-direction (while the structure in 
the Z-direction is clearly seen) , see Fig. [31 meaning that 
the core is loosing its crystalline structure locally and 
the Peierls potential in the climb direction is negligible 
under the simulated conditions. A systematic numeric 
study of the Peierls gap effects emerging at much lower 
temperatures and in larger system sizes remains a major 
computational challenge. 

Crucial data can be obtained experimentally with the 
"UMass sandwich" setup, that pot entially allows one to 
work at T of few tens of mK [25|. Since the quantity 
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FIG. 3: (Color online) A columnar snapshot of atomic posi- 
tions (open red dots) in the vicinity of the partial core (at 
the center) superimposed with the map (solid blue dots) of 
winding exchange cycles responsible for superfluid properties 
along the core, X-axis. Note that (i) the map extends over 
several unit cells (« 3.67A ), and (ii) it has no visible structure 
in the Y-direction, implying negligibly small Peierls potential 
for climb at the simulated temperature. 



of interest is the isochoric compressibility (as a function 
of T), one can use the superfluid syringe experimental 
protocol, when both Vycor electrodes are being operated 
at one and the same chemical potential and are used 
exclusively to inject atoms into the solid, rather than 



to induce a DC flow. Such measurements near 4007717^ 
have been already done [l^ . 

Summarizing, we present strong ab initio evidence and 
a coarse-grained analytic description of the climbing of an 
edge dislocation in solid ^He, assisted by superfluidity of 
its core. This phenomenon yields a natural microscopic 
interpretation for the effect of anomalous isochoric com- 
pressibility accompanying superflow in the experiment 
by Ray and Hallock. Theoretically, we argued that at 
low T, the superclimb, and, correspondingly, the effect of 
anomalous isochoric compressibility, must be suppressed 
due to a crossover to a smooth dislocation. Experimen- 
tal observation of the suppression, feasible within the 
"UMass sandwich" setup, might yield strong support for 
the proposed scenario bridging "quantum metallurgy" 
and supcrsolidity. The superclimb effect can also lead 
to high mobility of small dislocation loops (with Burgers 
vectors along C-axis) made of one partial surrounding 
an fee fault. Such loops could be plenty in real samples 
(cf. and imphcations of their presence are yet to 

be investigated. 
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